How occasional backstepping can speed up a processive motor protein 
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Fueled by the hydrolysis of ATP, the motor protein kinesin literally walks on two legs along the 
biopolymer microtubule. The number of accidental backsteps that kinesin takes appears to be much 
larger than what one would expect given the amount of free energy that ATP hydrolysis makes 
available. This is puzzling as more than a billion years of natural selection should have optimized 
the motor protein for its speed and efficiency. But more backstepping allows for the production 
of more entropy. Such entropy production will make free energy available. With this additional 
free energy, the catalytic cycle of the kinesin can be speeded up. We show how measured backstop 
percentages represent an optimum at which maximal net forward speed is achieved. 

PACS numbers: 05.40.-a, 87.16.Nn 

Processive motor proteins are among the tiniest engines known to man. These proteins utilize the energy of ATP 
hydrolysis to literally walk along a biopolymer pQ. In a living cell they help maintain organization by transporting 
cargo, like organelles or vesicles filled with chemicals. 

Already one and a half decade ago the stepping of the processive motor protein kinesin was made visible on the 
nanometer scale with optical tweezers pQ. Early communications [51 [3] reported that 5% to 10% of all steps of 
kinesin were backward. But smaller fractions were described later on as methods and materials improved and better 
resolutions were achieved; Ref. (4] gave 1/220 and Ref. [5] gave 1/802. Theoreticians have always been interested in 
backstep fractions as they can help verify stochastic models. However, throughout the literature backstepping has 
implicitly and explicitly been seen as an occasional malfunction of the stepping motor protein. 




FIG. 1: An abstract conception of kinesin's catalytic cycle. The cycle involves the binding and hydrolysis of ATP, the actual 
mechanical stepping, and the release of ADP and inorganic phosphate. Each reaction is driven by an energy difference Gi,i+i 
and the ratio of the forward and backward transition rate is an exponential function of G^i+i. 

In this Letter we will show how, in the Brownian environment of the motor protein, a "well-tuned" backstep fraction 
can actually help the motor speed up. We will show how the backstep fraction that leads to the highest net speed 
can be evaluated and how the resulting expression contains no freely adjustable parameters. Finally, we will see how 
the experimentally established backstep fraction of kinesin is close to our predicted optimal backstep fraction. 

The operation of an ion pump is generally modeled with a cycle as depicted in Fig. 1. At equilibrium the product 
of the forward rates, k\i x fc23 x ... x k n i, equals the product of the backward rates, foi x &32 x ... x k± n , and no 
net cycling occurs. To drive the protein through the sequence of states, Si, S2 ■■■ S n , a driving energy is necessary 
[5J. Such energy comes available if one of the steps involves the binding of ATP and if the protein, in subsequent 
steps, catalyzes the hydrolysis of the bound ATP. Eventually the remaining ADP and an inorganic phosphate have 
to be released so as to complete the cycle and to put the protein again in a state in which it can bind a new ATP. 
Under physiological conditions the hydrolysis of ATP makes Gatp = 22 fc^T units of free energy available. In the 
course of a cycle of a membrane pump like Na,K-ATPase, part of Gatp is utilized to bind, transport and release on 
the other side of the membrane three sodium ions and two potassium ions. The transport is generally against the 
electrochemical potential of the involved ions. Consistent with the model of Fig. 1, it is found that with a sufficiently 
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low ATP-ADP potential and a high electrochemical potential for sodium and potassium the operation of the pump 
can be reversed [7]- 



FIG. 2: The Brownian step of a processive motor protein. After detachment of the trailing leg, the attached leg reorients and 
brings the detached leg to the vicinity of the next forward binding site. After random diffusive motion will make the detached 
leg hit that binding site, actual binding can occur and a next step can commence. A backstep occurs after rebinding at the rear 
site that the detached leg just came from. The probabilities pj and pt for forward and backward binding depend on the energy 
behind the reorientation, the applied load force, and, as the process occurs in the Brownian regime, on ksT, the natural unit 
of thermal energy. 

The stepping kincsin appears to operate in a fundamentally different way. Tight coupling, i.e. an 8 nm step for 
every hydrolyzed ATP and a hydrolyzed ATP for every 8 nm step, has been observed for kinesin [5J [3] . Without a 
mechanical load it is just the Gatp that is driving the cycle in Fig. 1. Every mechanical step should correspond to one 
revolution around the catalytic cycle. If a backstep would correspond to the cycle being run in the direction against the 
ATP hydrolysis, then we should have Pb/pf = exp[— Gatp] for the ratio of the backstep probability and the forward 
step probability. However, exp[— 22] turns out 7 orders of magnitude smaller than the measured backstep fractions 
that were mentioned before. Furthermore, a model like in Fig. 1 leads to a stopping force, i.e. the load at which the 
kinesin comes to a standstill, that is determined by Gatp = FstL, where L is the stcplcngth. If Gatp = F st L, then 
the two forces, chemical and mechanical, that are driving the cycle in opposite directions cancel each other out. But, 
with Gatp = 22 and L = 8 nm, the equation Gatp = F st L predicts a stopping force F st that is about four times as 
large as the measured 7 pN [2l|3j[5]. Most importantly, it appears that kinesin still hydrolyzes ATP when it is pulled 
back with the stopping force and even when it is made to walk backwards with a load larger than the stopping force 
[5J[TU]. All these observations make a model as depicted in Fig. 1 untenable. 

We are thus led to a different model for the stepping motor protein [11]. After the detachment of the trailing leg, 
the attached leg reorients and brings the detached leg to the vicinity of the next forward binding site (see Fig. 2). 
Brownian motion is eventually supposed to make the detached leg hit the next forward site. Binding can then occur 
and a next step can commence. However, the energies involved in this Brownian step are comparable to fc^T. The 
energy G that drives the reorientation of the attached leg is smaller than Gatp- The detached leg can therefore 
rebind at the rear binding site that it came from and this can trigger the observed backstep. We can thus arrive at 
Pb/pf 3> exp[— Gatp] |12j . The scheme depicted in Fig. 3 is a more appropriate model for the Brownian stepper. At a 
particular point in the ATP hydrolysis cycle a kind of coin-toss takes place and the forward-backward determination 
occurs. The corresponding mechanical steps run in a dimension that is perpendicular to the plane of the chemical 
cycle. 

The question that needs answering at this point is: why has natural selection led to a backstep probability that 
is many orders of magnitude larger than the bare minimum of exp[— 22] that thermodynamics requires? As was 
mentioned before, kinesin pulls organelles or chemical-filled vesicles across a cukaryotic cell. The reason that eukaryotic 
cells have such an active transport system in the first place is that they are, unlike prokaryotic cells, too large to 
rely on diffusion for their transport needs. Kinesin 's speed ultimately determines how fast a eukaryotic cell can react 
to environmental stimuli. There is an obvious selectional advantage in engaging a kinesin that runs faster. No such 
straightforward evolutionary pressure for speedy turnover exists in the case of the aforementioned membrane pumps. 
A cell or organelle can meet its need to extrude sodium ions by simply putting the right amount of Na,K-pumps in 
the membrane. 

Our claim is that the entropy that is involved in the forward versus backward "choice" can provide an answer to 
the question. The free energy fc^T In 2 that is associated with a doubling of the available space is still small relative to 
the 22 fcsT-units released by ATP hydrolysis, but, we will argue, it can play a role in the optimization of a processive 
motor protein 
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FIG. 3: Kinesin still hydrolyzes ATP when it is pulled back with a force larger than the stopping force. A setup like Fig. 
1, with the ATP-ADP potential and the mechanical load pushing a single cycle in opposite directions, is therefore not the 
appropriate model and needs modification. Here we let the mechanical dimension run perpendicularly to the plane of the 
chemical cycle. ATP hydrolysis drives the chemical cycle in a clockwise direction. At the mechanical junction (cf. Fig. 2) a 
forward vs. backward "decision" is made. 

Imagine a large number, N, of motor proteins that are all at the same position. There is only one possible microstate 
for this macrostate. Next, let all of these N motors take a step. If all of these motors step forward, we will continue 
to have only one microstate for the macrostate of the system. But if, on the other hand, we allow for one backstep 
in the N steps, then there are N possible microstates. This is because the backstep can occur with any of the N 
motors. The increase in the number of microstates implies an increase in entropy. Another way to think of this is as 
follows. Before a step occurs, let the entropy be So = ks^a^lo, where f2o represents the number of microstates. After 
one catalytic cycle and its accompanying step the number of microstates regarding the positions of the motors on the 
microtubules is The internal states of the motors are the same after the cycle as before, so the total number of 
microstates after the step is f2i = f2o^p- For the ensuing entropy we have Si = ks In 17 = So + feeing. The total 
increase in entropy is due solely to the increase in "position entropy." We thus have AS = — ks (pfhipf + pblnpb) 
for the entropy increase that is associated with one step [13] . 

With a forward step probability p / and a backstep probabilty pb we have for the speed of the stepper v oc pf — pt = 
1 — 2pb- So the direct mechanical effect of an increase of the backstepping rate is a decrease of the speed, i.e. 
Sv oc — 2 Spb- But the backstep probability pb also makes free energy available to the amount of: 



Here G is expressed in units of ksT. ksT can be thought of as the natural unit of thermal energy and we will express 
energy in fe^T-units throughout the remainder of this Letter. As pb is observed to be small (pb < 0.01), first order 
approximations will be sufficiently accurate and the theory we develop in this Letter is a first order theory in pb- At 
first order in pb, Eq. M reduces to G w Ph(l — Inpb). So a variation 6pb in pb leads to a variation 



in the free energy that becomes available. 

We contend that the G of Eq. |l]) can be utilized to speed up the catalytic cycle of the motor protein. Throughout 
biology, energy is used to build and maintain ordered, low entropy structures. The sodium-glucose cotransporter 
concentrates glucose inside a cell. It does so by using the energy of a cotransported sodium ion. Entropic energy, 
in turn, is utilized whenever a transporter exploits a concentration gradient as a source of energy. In this way there 
are many examples of biomolecules that use the generation of entropy as a source of free energy. In the course of its 
cycle the Na,K,2Cl-cotransporter brings one sodium ion, one potassium ion, and two chloride ions from the outside to 
the inside of a cell. The cycle is clectroneutral. Under physiological conditions, this cotransporter effectively uses the 
transmembrane concentration gradients for sodium and chloride as an energy source to accumulate potassium inside 
the cell. 

We next need to quantify how a SG (cf. Eq. ^) can speed up a catalytic cycle. To this end we return to the 
catalytic cycle depicted in Fig. 1. The chemical part of kinesin's catalytic cycle (the binding of ATP, release of ADP 
and inorganic phosphate, detachment and attachment of the legs to the microtubule) is well modeled by a setup like 
Fig. 1. With energy expressed in units of fc^T, the energy that drives the transition from S, to Si+± is G^i+i and we 
have expfG^j+i] = fei j+i/A^+n PQ. Adding SG to G^i+i we can achieve an increase Ski^+i in the forward rate. 



G = TAS = -(pflnpf +p b lnp b ). 
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For small SG we have exp[5G] « 1 + SG, which leads to 

SG^ 6 ^. (4) 

So, ultimately, SG is just the relative increase of fc^i+i. It is also possible to put the SG towards reducing the reverse 
rate fcj+i^. However, when a reaction is already reasonably irreversible (i.e. fc»,i+i/fci+i,i being sufficiently large) this 
is not an efficient way to speed up the reaction. We thus assume that the SG brings down the energy of the product 
state and the energy of the activation barrier by the same amount. In the following we will ignore fti+i^ and write fcj 
for fci,i+i. 

The energies associated with the measured p b s are small. For p b = 1/802, Eq. (J!} gives G = 0.01 fcsT-units. 
Following Eq. Q we see that this translates into a 1% change in the speed of a transition relative to p b = 
case. This may seem insignificant. However, as we pointed out before, fast intracellular transport is an important, 
straightforward, and invariable selection criterion that is independent of the fitness landscape. 

When one transition in a cycle like in Fig. 1 is speeded up by 1%, it does not imply that the time to go through 
the entire cycle will decrease by 1%. That latter change will generally be less than 1%. We let C£. be a so-called 
"control coefficient" [HI [T3] ; this is a dimensionless parameter that puts a figure on the influence that the rate fe« of 
transition i has over the cycling rate v through the entire cycle. This cycling rate is also the stepping speed of the 
motor protein. We take 



Sv/v 
Ski/ k{ 



ci = ^r- (5) 



Going to the limit Ski ~* 0; we see that the control coefficient is essentially a logarithmic derivative, i.e. 
(dhiv) / (dlnki). So suppose fcj is changed by 1%. The control coefficient C£. then gives roughly the percentage 
by which v changes as a result. Assume that you change all of the transition rates in the cycle with the same percent- 
age. This is like scaling the time and, obviously, the speed v would be changed by the same percentage. This leads to 
the identity ^fc = 1' which is commonly known as the Summation Theorem |141 115] . If fcj represents the rate 

limiting transition, then C£. will be close to unity. C£. will be close to zero if the fc^-transition is very fast compared to 
other transitions in the cycle. Generally, the control coefficient C%. will be a positive number between zero and unity 
that expresses the amount of control that step i has over the speed of the entire cycle. In most biochemical networks 
and protein catalysis cycles there is no single rate limiting step. Instead, control is somewhat equally distributed over 
several steps. 

Since SG already denotes the relative change of fcj, we have for the speed of the motor protein and its increase due 
to the addition of SG to step i 

v + Svozl + Cl.SGttl-Cl. Qnp b )S Pb . (6) 

We combine this result with the mechanical effect of the backstepping, which, we saw earlier, was described by 
v oc 1 — 2p b or, consequently, 5v cx — 2 5p b . We now have for the net variation of v when p b is varied: 

Sv cx - {2 + CI lnp 6 } Sp b . (7) 

The speed v has a maximum when the variation Sv equals zero, i.e. 

lnp b = -2 /CI . (8) 

This equation relates the optimal backstepping fraction to the control coefficient of the transition that includes the 
free energy derived from the forward- versus-backward entropy. The p b = 1/220 from Rcf. [2] leads to Ct. — 0.4. The 
p b = 1/802 from Ref. [5] leads to C£ = 0.3. 

Under physiological conditions the transitions in the cycle of a motor protein are driven by energies of about 2 k^T- 
units |121 116] - Consequently, reverse transitions in the catalytic cycle are relatively rare. If we ignore the backward 
rates in the cycle of the motor protein, then the control coefficient CZ. is simply the fraction of time during the cycle 
that the motor protein spends in state Si. That time equals l/fej^+i. When the back leg of kinesin detaches, the 
hydrolysis of the ATP bound to that leg proceeds concurrently. The release of the ADP is supposed to take place 
upon rebinding of the detached leg to the microtubule. In a comprehensive review |16] a rate of ~ 250 s _1 is given for 
the back-leg-detachment /hydrolysis transition. With the stepping rate of ~ 100 s _1 , we are indeed led to a control 
coefficient of 0.4 for this transition. The same reference lists ~ 300 s _1 as the rate for the ADP release transition. 
This would correspond to a control coefficient of about 0.3. Both numbers are consistent with the backstepping rates 
of Refs. [3) and [5] and Eq. ^ and its underlying hypothesis. 
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There is no specific biomolecular mechanism that we can point at for the utilization of the free energy that is 
made available by a slightly increased pb and, in that sense, a solid underpinning of our hypothesis is lacking. It is 
generally true that having more configurations available within a macrostate effectively lowers the free energy of that 
state. An example illustrates how this is often important on the level of a single biomolecule. For a protein there are 
many ways to be unfolded and randomly coiled, i.e., there are many, say f2 u , unfolded microstates. The number of 
microstates, 0/, for the correctly folded state is generally orders of magnitude smaller [17 . So the folding of a protein 
actually requires a free energy TAS = kgTlii(fl u /flf). By making the folded state more flexible, i.e. by giving it 
more "statistical wiggle room," this required free energy may be brought down and that may lead to faster folding. 
To a large extent, the situation is analogous for the processive motor protein. By putting more available microstates 
at the end of the ATP hydrolysis cycle the motor can increase the free energy that is available to drive the step. 
Molecular dynamics simulations have been widely used to shed light on the dynamics of protein folding. They may 
become similarly illuminating for the case of stepping motor proteins. What we have shown in this Letter is that for 
kinesin the measured percentage of backsteps and the measured rates for the conformational changes agree with the 
idea of occasional backstepping as a source of energy and a way to speed up the motor protein. 
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